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' Abstract. Merging neutron star (NS) binaries may be detected by ground-based gravitational 

I ' wave (GW) interferometers (e.g. LIGO /VIRGO) within this decade and may also generate elec- 

, tromagnetic radiation detectable by wide-field, fast imaging telescopes that are coming online. 

^SJ ' The GWs can provide new constraint on the NS equation of state (including mass-radius relation 

^ ' and the related nuclear symmetry energy through resonant g-modes). This paper reviews vari- 

^ ous hydrodynamical (including equilibrium and dynamical/resonant tides) and electrodynamical 



0^ 



processes in coalescing NS binaries, with focus on the pre-merger phase. 
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\ 1. Introduction 

X : I was asked to talk about thermal radiation from isolated neutron stars (NSs). In this 

■ meeting, George Pavlov reviewed the X-ray properties of pulsars and thermally emitting 
NSs (see Kaplan et al. 2011), and Wynn Ho discussed central compact objects and their 

Q ■ magnetic fields (see Halpern & Gotthelf 2010; Shabaltas & Lai 2012; Vigano & Pons 
' 2012). Recent works on theoretical modelling of NS surface emission can be found in 

■ Potekhin et al. (2012) (see also Pavlov et al. 1995; Harding & Lai 2006 and van Adelsberg 
& Lai 2006 for reviews). Since these subjects were adequately covered in the meeting, I 
decided to focus on a different topic that did not receive much attention in this meeting 

T-H , but is likely to become increasingly important in the coming decade. 
1^ ■ Merging NS binaries have been studied since 1970s, with major activities in the rela- 
tivity community since the early 1990s because of their importance as a source of grav- 
itational waves (GWs) (e.g. Cutler et al. 1993). They are of great current interest for 
^ two reasons: (i) Merging NS/NS or NS/Black-Hole (BH) binaries have been identified 
as the leading candidate for the central engine of short GRBs (Berger 2011). They are 
I also expected to produce optical and radio transients that may be detected by wide-field, 
fast imaging telescopes that are coming online (e.g. PTE, LSST) in the next few years 
(Nissanke et al. 2012). (ii) After several decades of developmemt and promise, gravita- 
tional wave astronomy in the Hz-kHz band may finally take off in the next decade. The 
initial LIGO reached the design sensitivity [he — 10"^^) in 2006, and the enhanced LIGO 
/\ • (with a factor of 2 reduction in he) is taking or analysing data. The Advanced LIGO and 
VIRGO are expected to begin observations in 2015 and reach full sensitivity (a factor 
" " ■ of 10 reduction in he) in 2018-19 — at which time the detection of GWs from many 
merging NS binaries seems guaranteed. 

The last three minutes of a NS binary's life may be divided into two phases: the 
inspiral phase, producing quasi-periodic GWs, and the coalescence phase, where physical 
collision results in "messy" GWs. The recent years, 3D simulations of the final merger 
in full general relativity (GR) have become possible (see Shibata & Taniguchi 2006; 
Foucart et at. 2012; Sekiguchi et al. 2012). It has long been recognized that the final 
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merger waveforms can provide a useful probe of NS equation of state (EOS; e.g., Cutler 
et al. 1993; Bildsten & Cutler 1992; Lai & Wiseman 1996; Wiggins & Lai 2000). The 
idea is simple: By measuring the "cut-off" frequency oc {GMt/ R^Y^"^ associated with 
binary contact or tidal disruption, combined with the precise mass measurement from 
the inspiral waveform, one can obtain the NS radius (see Bauswein et al. 2012 for recent 
simulations which put such a idea into concrete footing; see also Sekiguchi et al. 2012; 
Faber & Rasio 2012 for reviews). 

In the following sections I will focus on the pre-merger phase. 

2. Hydrodynamics of merging NS binaries 

Prior to binary merger, tidal effects may affect the orbital decay and the GWs. There 
are two types of tides: equilibrium tides and dynamical tides. The equilibrium tides corre- 
spond to global deformation of the NS, which leads to the interaction potential between 
the two stars (with the NS mass M and radius R, the companion mass M' - treated as 
a point mass, and the binary separation a) 

V{r) = -MM' /a - O (kiM'^R^/a^^ , (2.1) 

where k2 is the so-called Love number. This would lead to a correction to the num- 
ber of GW cycles, dN = dN^°'> [l - 0{k2M'R^/Ma^)] . For a Newtonian polytropic NS 
model, simple analytic expressions can be found in Lai et al. (1994). Recent semi-analytic 
GR calculations of such equilibrium tidal effects (including the more precise determina- 
tion of the Love number) can be found in numerous papers (e.g., Flanagan & Hinderer 
2008; Binnington & Poisson 2009; Damour & Nagar 2009; Penner et al. 2012, Ferrari 
ct al. 2012). Obviously this effect is only important at small orbital separations (just 
prior to merger) - there is some prospect of measuring this, thereby constraining the 
EOS, but it may be challenging (Damour et al. 2012). At small orbital separations, the 
quadrupole approximation is not valid; therefore one must use the numerically computed 
GR quasi-equilibrium binary sequences to characterize the tidal effect such sequences 
have been constructed by several groups since the 1990s (e.g., Baumgarte et al. 1998; 
Uryu ct al. 2009) or use fully numerical simulations. 

Another aspect of the equilibrium tide concerns tidal dissipation, which leads to a lag 
of the tidal bulge with respect to the binary axis. It was shown already in the 1990s 
(Bildsten & Cutler 1992; Kochanek 1992) that because of the rapid GW-driven orbital 
decay, viscous tidal lag cannot synchronize the NS spin. Thus the NS will be close to 
irrotational (approximated as a Riemann-S ellipsoid; Lai et al. 1994; Wiggins & Lai 2000; 
Ferrari et al. 2012). Near the final phase of the inspiral, the rapid orbital decay gives rise 
to a finite lag angle (even with zero viscosity), but this cannot synchronize the NS (Lai 
& Shapiro 1995; Dall'Osso & Rossi 2012). 

The situation is more complicated for dynamical tides, which manifest as resonant 
excitations of internal oscillations of the NS: As two NSs spiral in, the orbit can momen- 
tarily come into resonance with the normal modes (frequency w^) of the NS: 

= mfiorb, m = 2, 3, • • • (2.2) 

By drawing energy from the orbital motion and resonantly exciting the modes, the rate 
of inspiral is modified, giving rise to a phase shift in the gravitational waveform. This 
problem was studied by Reisenegger & Goldreich (1994), Lai (1994) and Shibata (1994) 
in the case of non-rotating NSs, where the only modes that can be resonantly excited 
are g-modes (with typical mode frequencies;^ 100 Hz). It was found that the effect is 
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small for typical NS parameters (mass M = 1.4Mq and radius i? = 10 km) because the 
coupling between the g-mode and the tidal potential is weak. Ho & Lai (1999) studied 
the effect of NS rotation, and found that the g-mode resonance can be strongly enhanced 
even by a modest rotation (e.g., the phase shift in the waveform A$ reaches up to 
0.1 radian for a spin frequency ^ 100 Hz). They also found that for a rapidly rotating 
NS [vs <; 500 Hz), f-mode resonance becomes possible (since the inertial- frame f-mode 
frequency can be significantly reduced by rotation) and produces a large phase shift. In 
addition, NS rotation gives rise to r-mode resonance whose effect is appreciable only for 
very rapid (near breakup) rotations. Lai & Wu (2006) further studied resonant excitations 
of other inertial modes (of which r-mode is a member) and found similar effects. Flanagan 
& Racine (2006) studied the gravitomagnetic resonant excitation of r-modes and and 
found that the post-Newtonian effect is more important than the Newtonian tidal effect 
(and that the phase shift reaches 0.1 radian for Vg ~ 100 Hz). Tsang et al. (2012) 
examined crustal modes and found that the GW phase correction is small/modest and 
suggested that tidal resonance could shatter the NS crust, giving rise to the pre-cursor 
of short GRBs. Taken together, these studies suggest that for canonical NS parameters 
[R ~ 10 km, Vs ^ 100 Hz), tidal resonances have a small effect on the gravitational 
waveform during binary inspiral. However, it is important to keep in mind that the effect 
is a strong function of R (e.g., A$ oc i?^ for g-modes and cx R^-^ for inertial modes). A 
larger radius (i? ~ 15 km), appropriate for stiff EOS, would make the effect important. 
In the case of g-modes, the magnitude of the effect depends on the symmetry energy of 
nuclear matter and could be non-negligible (W. Newton & D. Lai 2013, in prep). 

3. Electrodynamics of merging NS binaries 

For magnetic NSs, magnetic interactions may play a role. If the binary is embedded in 
a vacuum, then the interaction potential is V{r) = —MM' /a — 0(/i/i'/a^) (where 
are the magnetic dipole moments of the two stars). It is easy to check that such magnetic 
interaction would lead to negligible effect on the GWs unless both NSs have superstrong 
fields (> 10^^ G) - this is unlikely (e.g., the double pulsars PSR J0737-3039 has 10^° G 
for pulsar A and 2 x 10^^ for pulsar B). 

Of course, as in the case of isolated pulsars, the circumbinary environment cannot 
be vacuum. The following discussion is based on Lai (2012). Consider a binary sys- 
tem consisting of a magnetic NS (the "primary", with mass M, radius R, spin $7^, and 
magnetic dipole moment /i) and a non-magnetic companion (mass Mc, radius Rc)- The 




T 

Figure 1. DC circuit (unipolar induction) model of magnetic interactions in binary systems a 

la Goldreich & Lynden-Bell (1969). 
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orbital angular frequency is il. The magnetic field strength at the surface of the primary 
is Bi, = n/R^. The whole binary system is embedded in a tenuous plasma (magneto- 
sphere). We assume for simplicity that S7, ilg and fi are all aligned. The motion of the 
non-magnetic companion relative to the magnetic field of the primary produces an EMF 
£ ~ 2Rc\E\, where E = Vi-d x B/c, with Vj-d = (il — fls)a^ and B = {—fj,/a^)z. This 
gives £ c± (2^i?c/ca^)Afi, where Afl = — fiJJ. The EMF drives a current along the 
magnetic field lines in the magnetosphere, connecting the primary and the companion 
through two flux tubes. The current in the circuit is given by I = £/{TZtot), where the 
total resistance of the circuit is 7?.tot = TZ + TZc + 27^mag, with TZ, TZc, TZmag the resis- 
tances of the magnetic star, the companion and the magnetosphere, respectively. These 
resistances depend on the properties of the binary components and the magnetosphere, 
and can vary widely for different types of systems. The energy dissipation rate of the 
system is then -Ediss — 2I^i?tot — 2£^/7^tot, where the factor of 2 accounts for both the 
upper and lower sides of the circuit. 

The total magnetic force (in the azimuthal direction) on the companion is ~ 
{2Rc){2TBz/c), with Bz = — /i/a^. Thus the torque acting on the binary's orbital angu- 
lar momentum is T = Jorb — {i/cja RcIBz — — {4: fiRc/ca^) {£ /TZtot)- The torque on the 
primary's spin is IVls = —T (where / is the moment of inertia). The orbital energy loss 
rate associated with T is then Eovh = TVl.. 

The equations above show that the binary interaction torque and energy dissipation 
associated with the DC circuit increase with decreasing total resistance Tttot- Is there 
a problem for the DC model when T^tot is too small? The answer is yes. The current 
in the circuit produces a toroidal magnetic field, which has the same magnitude but 
opposite direction above and below the equatorial plane. The toroidal field just above 
the companion star (in the upper flux tube) is ~ —{2'n:/c)Jr, where Jr — —41/ {t^Rc) 
is the (height-integrated) surface current. Thus the azimuthal twist of the flux tube is 

= -B^+/Bz == 16vrci/ic^'R-tot), where Vrei = aAD, = a{D, - fi^). Clearly, when 7?.tot 
is less than IQvrei/c^, the flux tube will be highly twisted. 

Goldreich & Lynden-Bell (1969) speculated that the DC circuit would break down 
when the twist is too large. (For the Jupiter-Io system parameters, the twist \(^\ <C 1.) 
Numerous works have since conflrmed that this is indeed the case. Theoretical studies 
and numerical simulations, usually carried out in the contexts of solar flares and accretion 
disks, have shown that as a flux tube is twisted beyond ^ 1, the magnetic pressure 
associated with B^ makes the flux tube expand outward and the magnetic flelds open up, 
allowing the system to reach a lower energy state (e.g., Aly 1985; Lynden-Bell & Boily 
1994; Lovelace et al. 1995; Uzdensky et al. 2002). Thus, a DC circuit with C0 ^ 1 cannot 
be realized: The flux tube will break up, disconnecting the linkage between the two binary 
components. A binary system with 7?.tot ^ 16wrei/c^ cannot establish a steady-state DC 
circuit. The electrodynamics is likely rather complex, only a quasi-cyclic circuit may be 
possible (Lai 2012; see Aly & Kuijpers 1990): (a) The magnetic field from the primary 
penetrates part of the companion, establishing magnetic linkage between the two stars; 
(b) The linked fields are twisted by differential rotation, generating toroidal field from 
the linked poloidal field; (c) As the toroidal magnetic field becomes comparable to the 
poloidal field, the fields inflate and the flux tube breaks, disrupting the magnetic linkage; 
(d) Reconnection between the inflated fleld lines relaxes the shear and restore the linkage. 
The whole cycle repeats. 

f If the magnetic field does have time to fully penetrate the companion star due to the rapid 
orbital decay, then the EMF will be reduced by a factor of order the ratio of the skin depth and 
the stellar radius. I thank Anatoly Spitkovsky for pointing this out to me at the lAU conference. 
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In general, we can use the dimensionless azimuthal twist to parameterize the mag- 
netic torque and energy dissipation rate: 

T=l-aRlB,B^+ = -a^^, ^diss = -TAO = ^AO^. (3.1) 

The maximum torque and dissipation are obtained by setting ~ 1. If the quasi-cychc 
circuit discussed in the last paragraph is established, we would expect to vary between 
and ~ 1. Note that in the above, T is negative since we are assuming > fig. 

Gravitational wave (GW) emission drives the orbital decay of the NS binary, with 
timescale ^Q-yv — = 0.012 (a/30 km)^ s, where we have adopted M = 1.4Mq and 

mass ratio q = M^/M = I. The magnetic torque tends to spin up the primary when 
r2>r2s. Spin-orbit synchronization is possible only if the synchronization time tsyn = 
IQ,/\T\ is less than icw at some orbital radii. With / = kMR^, we find 



^syn 



"^(l^)(i)--°V(^)"(3oLr3, (3.) 

where on the right we have adopted k = 0.4 and R = Rc = 10 km. Clearly, even with 
magnctar-likc field strength {B^^ ^ 10^^ G) and maximum efficiency (C^ ^ 1), spin-orbit 
synchronization cannot be achieved by magnetic torque. For the same reason, the effect 
of magnetic torque on the number of GW cycles during binary inspiral is small. 
The energy dissipation rate is 



'diss 



^ (v.^ABlR'^Rlc ^^44^ ( \^( a \-W2 _^ 



where on the right we have used v^ci — o,^ (for Qs ^ ^) and adopted canonical pa- 
rameters (M = Mc = 1.4M©, R = Rc = 10 km). The total energy dissipation per In a 
is 

^ = ^...*ow . 8.9 X (^)'(^)-'''erg. (3,4, 

Some fraction of this dissipation will emerge as electromagnetic radiation counterpart 
of binary inspiral. It is possible that this radiation is detectable at extragalactic dis- 
tance. But this will depend on the microphysics in the magnetosphcrc, including particle 
acceleration and radiation mechanism (e.g., Vietri 1996; Hansen & Lyutikov 2001). 

If one assumes that the magnetosphere resistance is given by the impedance of free 
space, 7?.„iag = 47r/c, then the corresponding twist is = 2t;rei/ (ttc), which satisfies our 
upper limit. We then have 

• /VrelVBlR^R^^C ^r.Adf V/ a \-7 , 

This is in agreement with the estimate of Lyutikov (2011). 

The situation is similar for NS/BH binaries. In the membrane paradigm (Thornc et 
al. 1986), a BH of mass Mh resembles a sphere of radius Rc = Rh = 2GMh/c^ (neglect- 
ing BH spin) and impedance TZh = 47r/c. Neglecting the resistances of the magnetosphere 
and the NS, the azimuthal twist of the fiux tube in the DC circuit is = Av^ci/ ('"'c), which 
satisfies our upper limit. The energy dissipation rate is (cf. Lyutikov 2011; Mc Williams 
& Levin 2011) 
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where we have assumed Mbh/M ^ 1. Again, it is uncertain whether this radiation is 
detectable for binaries at extragalactic distances. 
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